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The genetic basis and experimental evolution of inbreeding 
depression in Caenorhabditis elegans 

IM Chelo\ S Carvalho\ M Roque\ SR Prouk^ and H Teotonio^ 

Determining the genetic basis of inbreeding depression is innportant for understanding the role of selection in the evolution of 
nnixed breeding systenns. Here, we investigate how androdioecy (a breeding system characterized by partial selfing and 
outcrossing) and dioecy (characterized by obligatory outcrossing) influence the experimental evolution of inbreeding depression 
in Caenorhabditis eiegans. We derived inbred lines from ancestral and evolved populations and found that the dioecious 
lineages underwent more extinction than androdioecious lineages. For both breeding systems, however, there was selection 
during inbreeding because the diversity patterns of 337 single-nucleotide polymorphisms (SNPs) among surviving inbred lines 
deviated from neutral expectations. In parallel, we also followed the evolution of embryo to adult viability, which revealed 
similar starting levels of inbreeding depression in both breeding systems, but also outbreeding depression. Under androdioecy, 
diversity at a neutral subset of 134 SNPs correlated well with the viability trajectories, showing that the population genetic 
structure imposed by partial selfing affected the opportunity for different forms of selection. Our findings suggest that the 
interplay between the disruptions of coevolved sets of loci by outcrossing, the efficient purging of deleterious recessive alleles 
with selfing and overdominant selection with outcrossing can help explain mixed breeding systems. 
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INTRODUCTION 

Determining the genetics of inbreeding depression is important to 
understand the role of selection in the maintenance of selfing and 
outcrossing within populations (Goodwillie et al, 2005; Jarne and 
Auld, 2006). When inbreeding depression is due to overdominant 
loci, selfing is disfavored because individuals that self produce more 
homozygous progeny than if they were to outcross (Ziehe and 
Roberds, 1989; Charlesworth and Charlesworth, 1990; Uyenoyama 
and Waller, 1991a). In contrast, selfing is favoured when inbreeding 
depression is due to deleterious recessive alleles (Lande and Schemske, 
1985; Charlesworth et al, 1990; Uyenoyama and Waller, 1991b). 
Linkage disequilibrium and other non-random associations among 
deleterious recessives can however result in 'associative' overdomi- 
nance (Ohta and Kimura, 1970; Ohta, 1971; Palsson and Pamilo, 
1999), thus confounding the effects of overdominant loci on the 
evolution of selfing (Ziehe and Roberds, 1989; David, 1999; Bierne 
et al, 2000). 

Loci underlying differentiation in local environments might also 
influence the evolution of selfing, especially if they coevolved together 
within populations and their disruption leads to outbreeding depres- 
sion (Lynch, 1991; Charlesworth et al, 1997; Epinat and Lenormand, 
2009). Selfing in these circumstances might be favoured because it 
reduces effective segregation and recombination (Nordborg, 2000). 
However, disruption of coevolved sets of loci can also expose genetic 
variation at other loci that causes selection for either outcrossing or 
selfing, depending on whether the newly exposed variation is 
characterized by overdominance or partial dominance (Nordborg 
et al, 1996; Schierup and Christiansen, 1996; David, 1999). 



There is a poor empirical understanding about whether a balance 
between different forms of selection can explain mixed breeding 
systems (GoodwiUie et al, 2005; Escobar et al, 2008; Chelo and 
Teotonio, 2013). To address this problem, we used experimental 
evolution to study the genetic basis of inbreeding depression in large 
Caenorhabditis elegans populations. These populations were charac- 
terized by having either the wild-type male-hermaphrodite andro- 
dioecious breeding system (Maupas, 1900; Stewart and Phillips, 2002) 
or by having a male-female dioecious breeding system (Teotonio 
et al, 2012). We measured the evolutionary response in the prob- 
ability of survival with increased inbreeding levels, viability and the 
diversity of single-nucleotide polymorphisms (SNPs). Our findings 
suggest that a multistep process involving different forms of selection 
is responsible for mixed breeding systems. 

MATERIALS AND METHODS 

Experimental evolution 

The construction of the populations and experimental evolution design has 
been previously detailed (Teotonio et al, 2012). Briefly, the ancestral 
androdioecious population (termed EEV-Aq) resulted from a funnel pairwise 
cross among 16 wild isolates, while the ancestral dioecious population (EEV- 
Dq) was derived by the recurrent introgression of the fog-2(q71 ) allele into the 
Aq population for an extra 22 generations. The fog-2(q71 ) allele knocks out 
self-sperm (Schedl and Kimble, 1988), transforming hermaphrodites into 
functional females without apparent consequences for male reproductive 
success (Teotonio et al, 2012). 

Ancestral populations defined generation zero (GO) and there was threefold 
replication for experimental evolution (Ai_3 and Di_3), Populations were 
cultured alongside for 100 generations at constant 20 °C and 80% relative 
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humidity, under discrete 4-day non-overlapping life-cycles at census sizes of 
N= 10^ (Teotonio et al, 2012). Population samples were periodically stored 

-80 °C 

Inbreeding assays 

Population samples were revived from —80 °C stocks, each with > 10^ 
individuals, and cultured alongside for two generations under common 
environmental conditions. In the third generation, L3-L4 larval-staged 
(immature) individuals were sampled for the inbreeding assays. 

Inbred lines of the androdioecious populations were derived by selfing of 
hermaphrodites for 10 generations, from the Aq ancestral population, and from 
each of the three replicate populations at generations 30 and 100. Dioecious 
inbred lines were derived from GO, G30 and GlOO populations by brother-sister 
mating for 20 generations, to ensure similar final inbreeding coefficients to 
those of the lines derived from the androdioecious populations (see below). 

Lineages were maintained in 12-well culture plates, fiUed with 3.5 ml of 
NGM-lite media and 5|il of O/N cultures of Escherichia coli HT115. L3- or 
L4-staged individuals were passaged every 4-7 days to new plates. If 
reproduction or survival of a lineage failed after 7 days, individuals from the 
previous transfer, kept at 4 °C, were allowed to reproduce to higher densities 
and the protocol repeated in the following passage. Extinction of a lineage was 
scored at the generation where passage was unsuccessfial, after three such 
attempts. 

GO samples were included together with G30 and/or GlOO samples in four 
blocks, defined by the common calendar date of the beginning of the 
inbreeding protocol. After inbreeding, lines grew to exhaust available food 
and were fiarther cultured in 9-cm Petri dishes for two generations at high 
densities and frozen at —80 °C. For the ancestral populations over 120 lineages 
were inbred, for androdioecious evolved replicate populations 72 lineages were 
inbred and for dioecious evolved replicate populations 48 lineages were inbred 
(Supplementary Table 1). More derivations were initiated for androdioecious 
populations because there was little extinction. 

Survival analysis 

Lineage survival with expected inbreeding coefficients were calculated using a 
Kaplan-Meier estimator with right-censored data (Therneau and Grambsch, 
2000). Expected inbreeding coefficients were defined as: /f = 1 —^^Hq, where 
t is generation of selfing or sib mating where extinction was scored, X a 
limiting rate quantity set to 0.5 for selfing and to 0.809 for sib mating, and Hq 
the average number of heterozygous genotypes before any inbreeding was 
done, as calculated in each replicate population (Crow and Kimura, 1970). Hq 
was previously reported for GO, G30 and GlOO population samples, at 334 bi- 
aUelic SNPs measured in chromosomes IV and X (Chelo and Teotonio, 2013). 

We tested for the differences in the risk of lineage extinction with either 
between breeding systems at each generation or between generations within 
each breeding system. Cox proportional hazards models were employed to 
calculate different risks of extinction at each breeding system, using block as 
strata (Therneau and Grambsch, 2000). Ties were handled with the Efron 
approximation. The formulation using the survival package in the R statistical 
software (R Development Core Team, 2006; Therneau, 2012): coxph(Surv 
(ft.extinction, censored.status) ~ strata(block) + mating or generation). Like- 
lihood ratio tests with 1 d.f were used. 

Genotyping of inbred lines 

Frozen stocks of GO and GlOO inbred lines were thawed and cultured alongside 
for two generations at high densities. In the third generation, 20-30 L3- or L4- 
staged individuals were sampled from each of the inbred lines. 

Genomic DNA from pooled individuals was prepared with the ZyGEM 
prepGEM Insect kit following the manufacturer's protocol (ZyGEM Corpora- 
tion Ltd, Hamilton, New Zealand). A total of 337 bi-aUelic SNPs along 
chromosomes IV and X were chosen from the genome sequence of the N2, 
CB4856 and CB4858 wild isolates, as previously described in the study by 
Chelo and Teotonio (2013). Information about these SNPs is found in 
Supplementary Table 2. Genotypes were obtained by mass determination, 
after PGR amplification and aUele-specific extension using the iPlex Sequenom 
MALDITOF platform (Bradic et at, 2011). A total of 26 genotyping runs were 



done, each incorporating a maximum of 380 different inbred lines. In each 
run, 1-4 SNP plexes were used. 

Quality control was performed on data including the genotypes from the 
experimental evolution populations (Chelo and Teotonio, 2013). We first 
excluded SNPs with >80% missing data across all samples followed by 
removal of the inbred lines with > 50% of missing SNP genotypes. After this, 
SNPs with > 10% of missing data followed by inbred lines with > 10% of 
missing genotypes were removed. 

Physical positions among SNPs were defined according to the C. elegans 
genome release WS220 (December 2010). Genetic positions among SNPs were 
obtained by linear interpolation for the two chromosomes, using the fianction 
approx in R, each defined with genetic sizes of 50 cM (Rockman and Kruglyak, 
2009). Sex determination in C. elegans is chromosomal with hermaphrodites/ 
females XX and males X0 (Hodgkin, 1987). For chromosome IV, SNPs were at 
densities of 9.4/100 kb (3.3 SNP/cM) and for chromosome X at densities of 
9.8/100 kb (3.5 SNP/cM). The number of inbred lines genotyped can be found 
in Supplementary Table 1 and sample size details per SNP in Supplementary 
Table 2. 

SNP diversity after inbreeding 

Genetic diversity among inbred lines was estimated with the previously 
ascertained SNPs from the study of Chelo and Teotonio (2013). Average 
effective number of haplotypes was calculated across windows of 10 SNPs with 
step sizes of 1 SNP along the genetic distance at each chromosome: /Zg = 
1/J^pf^, with pi being the proportion of haplotype i among inbred lines (Crow 
and Kimura, 1970). Mean pairwise SNP linkage disequilibrium was also 
estimated in 10 SNP windows as the composite identity disequilibria among aU 
four genotypes. A, assuming that they were the product of the gametic 
probabilities: = /S^lpaqaPh^h'^ with p and q being the proportions of the most 
and least common allele, respectively, of SNPs a and h (Weir, 1996). 

To compute h^ and r^, SNPs were first phased into haplotypes using 
fastPHASE 1.2 (Scheet and Stephens, 2006). For each sample of inbred lines 20 
random starts of the EM algorithm were employed with 200 haplotypes taken 
from posterior distributions. The number of clusters for cross-validation was 
set to 10 and SNPs with posterior probabilities of <0.9 were considered 
missing data. Note that this protocol accounts for within-population genetic 
structure and thus the reconstructed inbred lines were diploids that could 
contain two different haplotypes. 

Expected neutral genetic diversity after inbreeding 

Monte-Carlo simulations of selfing for 10 generations or full sib mating for 20 
generations were performed in order to provide the neutral credible limits on 
the observed h^ and r^. A total of 1000 simulations were conducted per 
replicate population with resulting haplotypes being sampled in the same 
numbers as the inbred lines. Chromosome IV and chromosome X were 
analyzed separately. Details on the simulation algorithm can be found in the 
study by Chelo and Teotonio (2013). 

Each run started by randomly sampling phased diploids from the experi- 
mental replicates, in equal numbers as those of the starting inbred lineages. 
Recombination was simulated by exchanging consecutive sets of alleles between 
the two parental haplotypes (defined as vectors of SNP alleles and ordered as in 
Supplementary Table 2). We assumed complete crossover interference and map 
sizes of 50 cM. Crossover occurred randomly between any two consecutive 
SNPs according to the probability given by the genetic distances between them. 
For fertilization, and in the case of inbreeding by selfing, two independent 
gametes were joined to obtain the individual progeny. For brother-sister 
mating, two genotypes were chosen and kept separately at each generation. For 
the X chromosome in particular, male genotypes were defined by a single 
haplotype, to reflect their X-nuU constitution, by including an extra sampling 
step after recombination in females. 

Viability of ancestral inbred lines 

Twenty-five inbred lines from each of Aq or Dq populations were randomly 
revived from —80 °C stocks and were cultured alongside for two generations 
under the same environmental conditions. In the third generation, at day 4 of 
the life cycle, 100 embryos were collected to 6-cm Petri dish plates, incubated 
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for 4 days under standard conditions and the number of adult offspring scored 
after this period. Per inbred line, three assay plates were set up with 
manipulations and scoring randomized across breeding systems. 

Viability was the proportion of adult offspring at each plate out of 100 
embryos. Mixed effects ANOVA models were done to estimate the differentia- 
tion among breeding systems (Venables and Ripley, 2002). Random inbred 
lines were modeled within each breeding system and differences estimated by 
REML with the lme4 package in R: (lmer(viability~ mating + (l|line)). For 
significance, we assumed that the estimated effects followed Student's 
t distributions with 1 d.f. 



Viability in outbred populations 

Viability assays were carried out in three separate blocks, each including the GO 
ancestral populations and one same-numbered replicate population of each 
breeding system from GIO, G20, G41, G70 and GIOO. For each block, revived 
population samples were cultured alongside for two generations under 
common environmental conditions. On the third generation, we set up the 
assays as above for the GO inbred lines. Five replicate plates were prepared per 
population sample. 

Differences among breeding systems were modeled at GO by ANOVA with 
fixed blocks and fixed breeding systems using the stats package in 
R: (Im (viability ~ block + mating)). Viability trajectories were separately ana- 
lyzed at each breeding system by fitting mixed effects ANOVAs and estimating 
differences among generations while considering random block: 
(lmer(viability~generation+ (l|block)). For significance, we assumed that 
the REML estimated effects followed Student's t distributions with 1 d.f 

To illustrate the shape of the trajectories we fitted two-segment regressions 
of the mean viability among replicates onto generation using the stats and 
segmented packages in R (Muggeo, 2009): (lm(viability'~ generation); 
segmented.lm(obj = lm_viability, seg.Z= ~ generation, psi = G40)). To test 
for the significance of the estimated slopes before and after the break-point 
Student's t distributions were assumed with 1 d.f. 



Evolution of SNP diversity in outbred populations 

For each replicate experimental population we estimated the mean individual 
heterozygosity and the deviations in single and multi-locus genotype propor- 
tions, from those expected with random mating and infinite sizes (here termed 
genotype identity disequilibria). All of these metrics were obtained at GO, GIO, 
G30, G70 and GIOO from a subset of the 334 SNPs from the study by Chelo 
and Teotonio (2013), to encompass only those 134 SNPs that are located in the 
intergenic regions of chromosomes IV and X. They cover 1/3 of the genome 
and have densities of 1.1 SNP per cM in chromosome IV and 1.6 SNP per cM 
in chromosome X (Supplementary Table 2). We assume that these SNPs were 
neutral markers to the loci under putative selection during experimental 
evolution. For GO, average sample sizes were of ~90 genotypes while for 
remaining generations average sample sizes were of ~41 genotypes, at each of 
the SNPs (Supplementary Table 3). 

Individual heterozygosity {Hj) is the proportion of heterozygous SNPs 
across both chromosomes within each individual. Hi is interpreted as the 
inverse of the expected IBD or the inverse of the expected inbreeding 



coefficient of a randomly sampled individual; c.f. pp.62-68 (Crow and 
Kimura, 1970). Single-locus genotype disequilibria were measured as the 
fixation index Fjs^ I —{HJH), with Hg being the observed heterozygosity 
across SNPs and H being the expected heterozygosity under Hardy- Weinberg 
proportions; c.f. pp. 104-108 (Crow and Kimura, 1970). For pairwise genotype 
disequilibria, we calculated the average covariance in pairwise SNP hetero- 
zygosity {g2)y as computed with the RMES software (David et a/., 2007; Jame 
and David, 2008). 

To measure the extent of gametic disequilibria we calculated the 'back- 
ground' value of 1^. Polynomial functions were first fitted for aU the SNPs used 
in the study by Chelo and Teotonio (2013) against genetic distance, separately 
at each of the six regions in the chromosomes IV and X that are known to have 
fairly constant recombination rates, as defined in the study by Rockman and 
Kruglyak (2009). The genetic distance at which 5% of the initial decay was 
reached in each of the six regions was calculated and the average of the 
intergenic SNPs pairwise combinations above this distance taken as the 
background LD {bkgLD). bkgLD is inversely correlated to effective recombina- 
tion rates, as is a function of gamete combinations (see above) and because 
selection should not distort genotype disequilibria among loci separated by 
large genetic distances. 



Correlations of viability with SNP diversity 

For selection to maintain excess diversity, as previously measured in the 
experimental outbred populations (Chelo and Teotonio, 2013), the mean and/ 
or the variance among individual inbreeding coefficients and genotype identity 
disequilibria must be associated with fitness (David, 1999; Bierne et al, 2000; 
Navarro and Barton, 2002; Szulkin et al, 2010). To illustrate these associations, 
the Pearson product moment correlations of viability with Hp Fjs, g2 and 
bkgLD were calculated using the means of the three replicate populations 
within each breeding system at GIO, G30, G70 and GIOO. Viability data for 
G20 and G41 were averaged per replicate in order to pair it with the SNP data 
at G30. For significance testing, Pearson coefficients were transformed to 
Fisher's z- coefficients and across aU generations one-tailed t-tests with 3 d.f. 
were performed. Spearman coefficients gave similar results (analysis not 
shown). 



RESULTS 

Survival upon inbreeding 

The proportion of surviving lineages with inbreeding by selfing in 
androdioecious populations or with sib mating in dioecious popula- 
tions is shown in Figure 1 (see also Supplementary Figure 1). In the 
ancestral generation, the survival rate of dioecious lineages was ^0.6 
after inbreeding (inbreeding coefficients of/^^ 1), corresponding to a 
2.5-fold increase in the probability of extinction over that of the 
androdioecious populations (LK ratio test = 72.1, logio P= —7.23, 
n = 255). These differences continued to be apparent during experi- 
mental evolution. At G30, dioecious populations had a twofold higher 
chance of going extinct when inbred than androdioecious populations 
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Figure 1 Survival rates with inbreeding. Proportion of lineages surviving multiple generations of inbreeding by selfing (red) or sib mating (blue) in the 
androdioecious or dioecious populations, respectively, at different generations of experimental evolution (GO, G30 and GIOO). Error bars show 2 x s.d. See 
also Supplementary Figure 1. 
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(LK ratio test = 56.8, logio P= —13.3, n = 360), a difference that by 
GlOO was reduced to a 1.5-fold higher extinction risk (LK ratio 
test = 45.1, logio P= -10.7, n = 360). 

Analysis of lineage survival at each breeding system across the three 
periods indicated no evolution under androdioecy (LK ratio 
test = 2.1, P= 0.148, n = 562) and a marginal increase during evolu- 
tion under dioecy (LK ratio test = 3.5, P = 0.063, n = 413). 



SNP diversity after inbreeding 

At GO, SNP diversity among the inbred lines deviated from that 
expected with neutral processes during inbreeding (Figure 2). In both 
androdioecious and dioecious ancestors, inbred lines had higher 
haplotype diversity (/Zg) at chromosome IV than expected (Figure 2a). 
For chromosome X, however, there was no excess diversity 
(Figure 2b). 

By GlOO, selfing androdioecious hermaphrodites resulted in 
increased in both chromosomes relative to neutrality and when 
compared with the diversity of the outbred populations from which 
the inbred lines were derived (Figures 2a and b). Inbreeding the 
dioecious populations achieved higher than neutral expectations 
but the inbred lines showed reduced diversity relative to that of the 
outbred populations. Further, in contrast to androdioecy, under 
dioecy higher hg was only apparent in the autosome. 

Linkage disequilibrium among the inbred lines (r^) generally 
followed neutral expectations despite breeding system, chromosome 
and generation of experimental evolution (Figure 2c and d). However, 
and particularly for GlOO, all measured tended to cluster by the 
lower credible limits of neutrality. In fact, two out of three replicate 
androdioecious populations had significantly lower than expected 
with neutrality. 
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Figure 2 Genetic diversity after inbreeding. Effective haplotype number 
(a, b) and linkage disequibrium (c, d) for chromosomes IV (a, c) and 
chromosomes X (b, d). Circles show the diversity measured after inbreeding 
among the androdioecious lines (red) or dioecious lines (blue). Points 
indicate the diversity of the experimental populations before inbreeding. 
Error bars show the 95% credible probability obtained with 1000 neutral 
simulations of inbreeding. 



Evolution of viability 

Before inbreeding, the ancestral dioecious population was 8% less 
viable than the ancestral androdioecious population (Figure 3a; 
Fi 26 = 9.89, P = 0.004; block n.s.). Both selfing and sib mating led 
to 13% reductions in the mean values observed among inbred lines, 
with androdioecious lines being 8% more viable than dioecious lines 
it= -4.7, P<0.001). 

There was evidence for the experimental evolution of viability only 
under dioecy. For dioecy, regression analysis showed that the break- 
point is at G60 (60.44 ± 8.61 s.d.; adj. = 81%), and both before and 
after slopes are significant (pre-G60 t=3A9, P<0.001; post-G60 
t= -4.08, P<0.001). For androdioecy, the break-point is at G20 
(19.14 ± 17.53 SD) although the model has a very poor fit (adj. 
r2=10%) and neither slope is significant (pre-G20 t= —0.728, 
P = 0.466; post-G20 t = 0.910, P = 0.363). 

Correlated evolution of viability with SNP diversity 

The correlations of SNP diversity with viability were most obvious 
under androdioecy (Figure 4). In this breeding system, tests across all 
generations showed a significant correlation of viability with bkgLD 
(P=0.02) and marginal significance with g2 (P=0.06). When not 
including GlOO, all androdioecious correlations are different from 
zero (P = 0.01). In contrast to androdioecy, under dioecy, tests across 
aU generations only revealed significance of the correlation of g2 with 
viability (P = 0.04). 

DISCUSSION 

Inbreeding and outbreeding depression 

Inbreeding depression is thought to mostly occur because deleterious 
recessive alleles are expressed in homozygotes (Charlesworth and 
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Figure 3 Evolution of viability, (a) Outbred or inbred viability of the 
ancestral androdioecious population (red) and the ancestral dioecious 
population (blue). Error bars show 1 s.e.m. among assay blocks for the 
experimental populations or one s.e.m. among the inbred lines, (b) Viability 
of androdioecious (red) and of dioecious (blue) populations during 
experimental evolution. Error bars indicate 1 s.e.m. among the three 
replicates. Lines show the estimated two-segment regressions. 
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Figure 4 Population structure and selection. Correlation coefficients of 
viability with several genotype disequilibria metrics during experimental 
evolution. In solid lines, comparisons that are significantly different from 
zero, gray dashed lines otherwise. 
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Willis, 2009). As selection against recessive deleterious alleles is weaker 
under outcrossing than under selfing, we expected that inbreeding 
depression would be maintained at higher levels under dioecy than 
under androdioecy (Lande and Schemske, 1985; Charlesworth et al, 
1990; Uyenoyama and WaUer, 1991b). As expected, we found that 
dioecious populations subjected to inbreeding had low lineage 
survival at all generations of experimental evolution, whereas andro- 
dioecious populations subjected to inbreeding had high lineage 
survival rates. 

The viability results in the ancestral populations are also consistent 
with inbreeding depression being due to deleterious recessives because 
the inbred lines were less viable than the outbred populations from 
which they were derived. Further, the viability data showed initial 
outbreeding depression, with dioecious populations having lower 
viability than androdioecious populations, regardless of the level of 
inbreeding. Interestingly, therefore, outbreeding and inbreeding 
depression were not independent phenomena; c./., (Lynch, 1991; 
Schierup and Christiansen, 1996; Escobar et al, 2008; Epinat and 
Lenormand, 2009). 

Outbreeding depression implies underdominance (Lynch, 1991), 
following the disruption by segregation and recombination of sets of 
loci that coevolved in different populations through adaptation to 
local environmental conditions or by genetic drift (Coyne and Orr, 
1998; Kirkpatrick and Barton, 2006; Epinat and Lenormand, 2009). 
However, underdominance would not have led to the excess diversity 
that we measured among the inbred lines. Instead, there must have 
been fitness overdominance during inbreeding. Notably, a similar 
conclusion was reached when excess diversity relative to neutral 
expectations was detected in the experimental outbred populations 
(Chelo and Teotonio, 2013). 

It is known that fitness overdominance can result from the non- 
random association of deleterious recessive alleles that are not 
necessarily in close physical linkage (Ohta and Kimura, 1970; Ohta, 
1971; Palsson and Pamilo, 1999). Could this associative overdomi- 
nance help explain the interdependence of outbreeding and inbreed- 
ing depression during experimental evolution? Specifically, as long as 
the coevolved sets of loci remained intact there would be outbreeding 
depression but once disrupted would there be inbreeding depression? 
An answer to this question needs a better understanding of how the 
population genetic structure imposed by the two breeding systems 
influences the opportunity for different forms of selection (David, 
1999; Szulkin et al, 2010). 

Population structure and selection 

As expected with littie population genetic structure (Charlesworth 
et al, 1990; Szulkin et al, 2010), under dioecy there was no trend in 
the correlations of SNP diversity with viability over the course of the 
experimental evolution, implying that selection was not very efficient 
at removing deleterious recessive alleles and/or at sustaining fitness 
overdominance. 

In the androdioecious populations, however, several correlations 
were significant. There was a negative correlation between viability 
with background linkage disequilibrium (bkgLD), which suggests that 
selection favoured new recombinants. This is because bkgLD should 
quantif)^ the extent of gametic linkage disequilibrium and thus 
effective recombination (Christiansen, 1989). Additionally, the corre- 
lations between the co-variation in diversity within SNPs (Fjs) and 
the CO -variation in diversity among SNPs {g2) with viability were 
positive. Positive signs in these correlations indicate that identity 
disequilibria at multiple loci across the genome facilitated selection 
among genetically heterogeneous individuals (David, 1999; Szulkin 



et al, 2010). This population genetic structure may have in turn 
reinforced selection against deleterious recessive alleles, which is 
consistent with the observation of a negative correlation of viability 
with individual heterozygosity (Hi). 

If the population genetic structure imposed by partial selfing in 
androdioecy enabled the purging of deleterious recessive alleles, then 
the fitness overdominance responsible for excess diversity in the 
outbred populations might have been due to truly overdominant loci 
(Christiansen, 1989; Ziehe and Roberds, 1989; Charlesworth and 
Charlesworth, 1990). In fact, the positive correlations between 
viability and Fjs and g2 found under androdioecy suggest selection 
on overdominant loci. This is because these metrics also quantify the 
number of heterozygote classes (Weir et al, 1980; David, 1999; David 
et al, 2007), and it has been theoretically shown that selection on a 
few overdominant loci creates positive correlations between the 
number of heterozygote classes and fitness variance (David, 1999). 
Furthermore, if the average population fitness was a diminishing 
returns function of heterozygosity the negative correlations of viability 
with bkgLD under androdioecy could be generated by overdominant 
loci (Navarro and Barton, 2002). 

Reductions in viability with inbreeding were accompanied by 
higher diversity among inbred lines than among the individuals of 
the outbred populations from which they derived, a pattern particu- 
larly evident in the autosomes at the beginning of experimental 
evolution. These results can be explained by the generation of new 
deleterious alleles during inbreeding, of which only some were 
selected against, c.f. (Barriere et al, 2009). Diversity results among 
inbred lines at generation 100 further support a role for truly 
overdominant loci. Dioecious sib mating led to less diversity among 
inbred lines than selfing under androdioecy. Despite the possibility for 
similar kinds of selection in the two breeding systems, upon 
inbreeding, reduced effective recombination with selfing could have 
resulted in higher (haplotype) diversity because of selection on 
overdominant loci (Navarro and Barton, 2002). 

Selection and the maintenance of androdioecy 

When sexual selection is considered the diversity results can be more 
fully explained (Anthes et al, 2010; Baer et al, 2010; Mallet and 
Chippindale, 2011). In particular, sexual selection should have been 
stronger under dioecy because higher numbers of males in this 
breeding system would have led to lower numbers of 
X chromosomes — in C. elegans males are X0 and hermaphrodites/ 
females are XX (Hodgkin, 1987). Following inbreeding there was 
excess diversity in the X-chromosome under androdioecy but not 
under dioecy. As expected, therefore, the sex ratio of a population 
might have been associated with the removal of deleterious recessive 
alleles from the X-chromosome. This result is remarkable because we 
were previously unable to confirm that sex ratio differences among 
breeding systems influenced the evolution of male competitive 
performance (Teotonio et al, 2012), presumably a fitness component 
under strong sexual selection, compare with (LaMunyon and Ward, 
2002; Murray et al, 2011). 

Taken together, our findings point to a multistep process that 
maintains partial selfing under androdioecy (also refer to 
Charlesworth et al, 1997; Pannell, 2002; Goodwillie et al, 2005). 
Transitions from outcrossing to selfing can first enable the appearance 
of different sets of coevolved loci, given sufficient time for differentia- 
tion among populations. These transitions to selfing could occur, for 
example, because of reproductive assurance during the colonization of 
empty habitats (Cheptou, 2004). However, some degree of out- 
crossing is inevitable because of recurrent mutation in sex 
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determination pathways and dispersal of males among neighboring 
populations. The disruption of coevolved sets of loci would in turn 
expose to selection partially dominant loci that originated during 
population differentiation, which would favor selfing. Because of 
purging deleterious recessives, however, transient associative fitness 
overdominance and/or selection on truly overdominant loci would 
allow outcrossing to persist until local population extinction and 
novel transitions to selfing. 

Several lines of evidence suggest the occurrence of this multistep 
process in the maintenance of androdioecy in C. elegans. Hybridiza- 
tion of wild isolates results in outbreeding depression that might be 
due to several loci in complete linkage disequilibrium (Dolgin et al, 
2007; Seidel et al, 2008). There is also abundant genetic variation for 
male function and selection for outcrossing in novel environments 
(Teotonio et al, 2006; Manoel et al, 2007; Murray et al, 2011; 
Teotonio et al, 2012), even if males are rarely found in natural 
populations (Felix and Duveau, 2012). Finally, heterozygosity within 
natural populations might be higher than that expected with a long 
history of exclusive selfing (but see Barriere and Felix, 2005; 
Sivasundar and Hey, 2005; Cutter et al, 2008; Andersen et al, 
2012). All these observations indicate that in C. elegans a balance 
between different forms of selection maintains both selfing and 
outcrossing. Perhaps a similar balance also explains mixed breeding 
systems in other species. 
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